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EFFICIENTLY DETERMINING A FLOOR FOR A 
FLOATING-POINT NUMBER 

BACKGROUND 

The present invention relates generally to computation using binary floating-point 
numbers^ and particularly to finding the floor of a floating-point number. 

The floor for a floating-point number x is the largest integer less than or equal to x So 
while the floor for 2.5 is 2, the floor for -2.5 is -3. Consequently, finding a floor is more 
difficult for negative floating-point numbers than for positive floating-point numbers. 

SUMMARY 

Li general, in one aspect, the invention features an apparatus, method and computer 
program product for processing a binary floating-point number having a sign bit and a 
mantissa having a fraction portion. It includes identifying the fraction portion of the binary 
floating-point number; and replacing each bit of the fraction portion with the sign bit, thereby 
producing a floor of the binary floating-point number. 

Particular implementations can include one or more of the following features. 
Implementations include decrementing the binary floating-point number before replacing 
when the binary floating-point number is negative. Implementations include converting the 
floor to two's complement format. Converting includes performing an exclusive-OR 
operation between each bit of the floor and the sign bit, thereby producing a result of the 
exclusive-OR operation; and concatenating the sign bit and the result of the exclusive-OR 
operation, thereby producing a signed two's complement mantissa of the floor. Converting 
further includes performing, upon the signed two's complement mantissa of the floor, a 
signed-right-shifl operation, thereby producing the floor of the binary floating-point number 
in two's complement format. 

Implementations include converting the floor to floating-point format. Converting 
includes incrementing the floor when the binary floating-point number is negative, and doing 
nothing otherwise, thereby producing an incremented value; and replacing the most- 
significant bit (MSB) of the incremented value with the exponent bits and the sign bit, such 
that the sign bit is the MSB, thereby producing a floor of the binary floating-point number in 
floating-point format. 
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The binary floating-point number includes an exponent that differs from an unbiased 
exponent by a bias offset, and incrementing includes incrementing the floor when the binary 
floating-point nxmiber is negative and the exponent is greater than, or equal to, the bias offset, 
thereby producing an incremented value. The binary floating-point number includes an 
exponent that differs from an unbiased exponent by a bias offset, and incrementing further 
includes replacing the exponent bits with the offset when the binary floating-point number is 
negative and the exponent is less than the offset; and replacing the exponent bits with zeros 
when the binary floating-point number is positive and the exponent is less than the offset. 

Implementations include taking a floating-point difference between a value of the 
binary floating-point number before replacing and a value of the binary floating-point 
number after replacing, thereby producing a fractional remainder of the binary floating-point 
number. 

The binary floating-point number includes an exponent that differs from an unbiased 
exponent by a bias offset, and replacing includes replacing each bit of the fraction portion 
with zero ("0") when the exponent is greater than, or equal to, the bias offset. The binary 
floating-point number includes an exponent that differs from an unbiased exponent by a bias 
offset, and decrementing includes decrementing the binary floating-point number before 
replacing when the binary floating-point number is negative unless the exponent is less than 
the bias offset. 

In general, in one aspect, the invention features an apparatus, method and computer 
program product for determining the floating-point floor of a floating-point number. It 
includes identifying a binary floating-point number including a sign bit, exponent bits, and 
mantissa bits, wherein the binary floating-point number is negative when the sign bit is a one 
("1"); concatenating an implicit bit and the mantissa bits, thereby producing a first binary 
number such that the implicit bit is the most significant bit (MSB) of the first binary number; 
decrementing the first binary number when the sign bit is a one (" 1 ") and doing nothing when 
the sign bit is a zero ("0"), thereby producing a second binary number; identifying a fraction 
portion of the second binary number based upon a predetermined exponent bias; replacing 
each bit of the fraction portion with the sign bit thereby producing a third binary number; 
performing an exclusive-OR operation between each bit of the third binary number and the 
sign bit thereby producing a fourth binary number; concatenating the sign bit, the fourth 



binary number, and a first predetermined number of zeros, thereby producing a fifth binary 
number such that the sign bit is the MSB of the fifth binary number, and the zeros are the 
least significant bits of the fifth binary number; and performing, upon the fifth binary 
number, a signed-right-shift operation by a number of bits equivalent to the difference 

5 between the exponent and a second predetermined number, thereby producing the floor of the 
binary floating-point number in integer format. 

Particular implementations can include one or more of the following features. 
Implementations include subtracting one from the simi of the number of bits in the floating- 
point floor and the number of bits in the third binary number, thereby producing the first 

10 predetermined number. The exponent differs fi*om an imbiased exponent of the floating-point 
number by a bias offset, and implementations include summing the bias offset, the number of 
bits in the mantissa, and the first predetermined number, thereby producing the second 

ri 

predetermined number. Replacing includes identifying a fraction mask corresponding to the 
^ exponent bits, the fi-action mask having a one for each bit belonging to the fraction portion 

Ifi 15 and a zero for each bit belonging to the integer portion; and applying the fraction mask to the 
n second binary number. The impUcit bit is zero ("0") when the exponent bits are all zero (**0") 

W and the implicit bit is one (" 1 ") otherwise. 

13 The binary floating-point number includes an exponent that differs from an unbiased 

^; exponent of the floating-point number by a bias offset, and replacing includes replacing each 

19 20 bit of the fraction portion with the sign bit when the exponent is greater than, or equal to, the 
bias offset, thereby producing a third binary number. The binary floating-point number 
includes an exponent that differs from an unbiased exponent of the floating-point number by 
a bias offset, and decrementing includes decrementing the first binary number when the sign 
bit is a one ("1") and doing nothing when the sign bit is a zero ("0"), unless the exponent is 
25 less than the bias offset, thereby producing a second binary number. 

Implementations include replacing the MSB of the second binary number with the 
exponent bits and the sign bit, thereby producing a sixth binary number, such that the sign bit 
is the MSB of the sixth binary number; replacing the MSB of the third binary number with 
the exponent bits and the sign bit, thereby producing a seventh binary number, such that the 
30 sign bit is the MSB of the seventh binary number; and performing a floating point subtraction 
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with the sixth binary number as the minuend and the seventh binary number as the 
subtrahend, thereby producing a fractional remainder of the floating-point number. 

Advantages that can be seen in implementations of the invention include one or more 
of the following. Implementations of the present invention determine one or both of the floor 

[x J and fractional remainder r = x- [x J for a floating-point number x. hnplementations that 

determine both the floor and the fractional remainder share common components, thus 
reducing the total number of components required for both functions. Further, 
implementations of the present invention employ a floating-point subtractor; because most 
general-purpose processors include a floating-point subtractor, such implementations can be 
incorporated in such processors with few additional components. 

The details of one or more embodiments of the invention are set forth in the accompa- 
nying drawings and the description below. Other features, objects, and advantages of the 
invention will be apparent from the description and drawings, and from the claims. 

DESCRIPTION OF DRAWINGS 

FIG. 1 is a logical block diagram of an apparatus according to an implementation of 
the present invention. 

FIG. 2 is a flowchart depicting a process executed by the apparatus of FIG. 1 to 
determine a floor for a floating-point number according to one implementation. 

FIG. 3 is a flowchart depicting a process executed by the apparatus of FIG. 1 to 
determine a fractional remainder for a floating-point number according to one 
implementation. 

FIG. 4 is a logical block diagram of a format converter according to an 
implementation of the present invention. 

FIG. 5 is a flowchart depicting a process executed by the format converter of FIG. 4 
to convert a floor from conditionally-decremented floating-point format to two's complement 
format according to one implementation. 

FIG. 6 is a logical block diagram of a format converter according to an 
implementation of the present invention. 



FIG. 7 is a flowchart depicting a process executed by the format converter of FIG. 6 
to convert a floor from conditionally-decremented floating-point format to floating-point 
format according to one implementation. 

Like reference symbols in the various drawings indicate like elements. 



described below for floating-point numbers defined by IEEE (Institute of Electrical and 
Electronics Engineers) Standard 754 for single-precision binary floating-point computations 
10 (hereinafter referred to as "IEEE 754"). According to IEEE 754, a floating-point number is 
described by a 32-bit binary number having one bit of sign s, eight bits of exponent e, and 23 
bits of mantissa m. A floating-point number is positive when ^ = 0 and negative when ^ = 1 . 
The value of the mantissa m lies between one and two such that 
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DETAILED DESCRIPTION 



Implementations of the present invention determine one or both of the floor \_xj and 
fractional remainder r = x - [xj for a floating-point number x One implementation is 



l<l + 2-'m<2 



(1) 



Therefore the floating point values ofx are given by 




x = i 



e = 0,m — 0 
e = 255,m = 0 



(2) 



GO 
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In one implementation, the input x is a 32-bit single-precision floating-point number 
according to IEEE 754, [x J is a 32-bit two's-complement integer, and r is a 32-bit single- 
precision floating-point number according to IEEE 754. However, the techniques of the 
present invention are applicable to inputs and outputs of different formats. 
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Implementations of the invention employ a biased exponent e that is simply a biased 
version of the exponent e' defined by IEEE 754. In the general case 
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(3) 



However, EEEE 754 defines special cases for certain values of exponent e. For e 0, IEEE 
754 specifies ^ -126, p ^ 0, and = 0. For the case e = 255, m = 0, and 5 = 0, IEEE 754 
5 specifies x= +co. For the case e = 255, m = 0, and ^ = 1, IEEE 754 specifies x = -00 . For the 
case e = 255 and m 1^ O.xis not a number (NAN); this enables special meanings to be 
assigned to this case. For the convenience of the reader, these relationships are shown in 
Table 1. 
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Table 1 also shows the range of the floating-point values of x in terms of the 
maximum {xmw) and minimum {xmitd floating-point values of x possible for each exponent e. 
Table 1 also depicts the implicit bit p for each exponent e. Note that p = 0 only when e = 0, 

5 and = 1 otherwise. 

Table 1 includes a fraction mask column for use with the present invention. In each 
row, Table 1 lists a value of exponent e and a corresponding fraction mask value c. For 
clarity. Table 1 depicts only the rows corresponding toe = 0, 125<e<131, 148<e<151, 
and e > 253. The contents of the remaining rows of Table 1 will be apparent to one skilled in 

1 0 the relevant arts after reading this description. 

The mantissa m includes one or both of a fraction portion and an integer portion. Thus 
each bit of the mantissa m may be part of the integer portion or the fraction portion. The 
fraction mask indicates which bits of the mantissa belong to the fraction portion. If a mask bit 
is 1, the corresponding mantissa bit is part of the fraction portion. If a mask bit is 0, the 

15 corresponding mantissa bit is part of the integer portion. Referring to Table 1, when e < 127, 
all of the mantissa bits belong to the fraction portion; when 126 < e < 150, some of the 
mantissa bits belong to the fraction portion and some to the integer portion; and when e > 
149, all of the mantissa bits belong to the integer portion. Note that the fraction mask is also 
applied to the implicit bit p. 

20 FIG. 1 is a logical block diagram of an apparatus 100 according to an implementation 

of the present invention. FIG. 2 is a flowchart depicting a process 200 executed by apparatus 
100 to determine a floor for a floating-point number according to one implementation. 

Process 200 identifies a binary floating-point number x including a sign bit s, 
exponent bits e, and mantissa bits m (step 202). In one implementation, floating-point 

25 number x is a 32-bit single-precision floating-point number x 

Process 200 removes from x the sign bit s and the exponent bits e, leaving only the 
23 -bit mantissa m (step 204). Process 200 then concatenates an implicit bit p and mantissa m 
(step 206), thereby producing binary number mj such that implicit bit p is the most 
significant bit (MSB) of binary number m;. Implicit bit p is selected according to the 

30 exponent e, as shown in Table 1. Such removing and concatenating of bits is accomplished 
by conventional wiring and bus routing structures apparent to one skilled in the relevant arts. 
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Process 200 decrements binary nximber m; by one if ^ = 1 (indicating that the 
floating-point number is negative), thereby producing binary number m2 (step 208); 
otherwise m2 = mj.ii the implementation of FIG. 1, the decrement operation is performed by 
a decrementer 102. In another implementation, as shown in FIG. 1, process 200 instead 
decrements binary number /w; by one if a conditional sign bit ^5 = 1, thereby producing 
binary number (step 208), where conditional sign bit ss is given by 

ife> 126,then55 = ^ (4) 
else S3 = 0 

Otherwise m2'^m}. 

Process 200 identifies a fraction portion of binary number m2 based upon the 
exponent e (step 210). In the implementation of FIG. 1, process 200 indexes a fraction mask 
table 106 using exponent e to obtain a fraction mask c. The fraction mask c has a one ("1") 
for each bit belonging to the fraction portion and a zero ("0") for each bit belonging to the 
integer portion. The fraction mask value c for each value of exponent e is given by Table 1 . 

Process 200 replaces each bit of the fraction portion with sign bit s. In another 
implementation, as shown in FIG. 1, process 200 replaces each bit of the fraction portion 
with conditional sign bit 55, thereby producing a third binary number ms (step 212). In the 
implementation of FIG. 1, process 200 appUes fraction mask c to binary number m2, one bit 
at a time, LSB first, thereby producing binary number m^. Each bit of fraction mask c is 
applied to a multiplexer 104 as the control input. The "1" input of multiplexer 104 is sign bit 
s. In another implementation, as shown in FIG. 1, the "1" input of multiplexer 104 is 
conditional sign bit 5-5. The bits of binary number m2 are fed sequentially to the "0" input of 
multiplexer 104. When a bit of fraction mask c is zero, multiplexer 104 gates the 
corresponding bit of binary number /w^- When a bit of fraction mask c is one, multiplexer 104 
gates sign bit s. In another implementation, as shown in FIG. 1, when a bit of fraction mask c 
is one, multiplexer 104 gates conditional sign bit ss. Binary number ms is a floor of floating 
point number ;c in a format referred to herein as "conditionally-decremented floating-point 
format." It is often convenient to convert floating-point floor ms to a more conventional 
format using a format converter 108 (step 214). 
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FIG. 3 is a flowchart depicting a process 300 executed by apparatus 100 to determine 
a fractional remainder r = jc - [x J for a floating-point number x according to one 

implementation. Process 300 replaces the MSB of binary number m2 with the exponent bits e 
and the sign bit s, thereby producing a binary number ms, such that the sign bit s is the MSB 
of binary number me (step 302). Similarly, process 300 replaces the MSB of binary number 
ms with exponent bits e and the sign bit s, thereby producing a binary number m?, such that 
the sign bit s is the MSB of binary number my (step 304). In another implementation, as 
shown in FIG. 1, process 300 replaces the MSB of binary number ms with conditional 
exponent bits ej and the sign bit s, thereby producing a binary number mj, such that the sign 
bit s is the MSB of binary number m?. Conditional exponent bits es are given by 



if e> 126 then 03 = 6 
else if 5- = 1 then e = 127 
else e^O 



(5) 



Process 300 then takes the floating-point difference m, -m^^ (step 306): that is, using 
binary number ms as the minuend and binary number my as the subtrahend. In the 
implementation of FIG. 1, the difference is computed by a floating-point subtractor 112. The 
difference is the fractional remainder r presented as a 32-bit single-precision floating-point 
number according to IEEE 754. 

FIG. 4 is a logical block diagram of a format converter 400 according to an 
implementation of the present invention. FIG. 5 is a flowchart depicting a process 500 
executed by format converter 400 to convert floor ms from conditionally-decremented 
floating-point format to two's complement format according to one implementation. 

Process 500 performs an exclusive-OR (XOR) operation between each bit of binary 
number ms and sign bit s, thereby producing 24-bit binary number m4 (step 502). In one 
implementation, the XOR operation is performed using an XOR gate 402 or some other such 
complementer. 

Process 500 concatenates binary number m4 and a predetermined number of zeros, 
such that the zeros are the least significant bits (LSB) of the result (step 504). The 
predetermined number of zeros z is given by 



z = bits(\_x\) - hits{m^ ) - 1 



(6) 



where ft/r.s(LxJ) is the number of bits in the floating-point floor \ x\ and hits{m^) is the 
number of bits in binary number mj. In the implementation of Table 1, hits(\_xj) = 32 and 

bits(m3) = 24, yielding z 7 . 

Process 500 concatenates the result and sign bit s, thereby producing a binary number 
ms such that the sign bit s is the MSB of binary number ms (step 504). Binary number mj is a 
signed two's complement mantissa of the floor of x. 

Process 500 performs, upon binary number mj, a signed-right-shift operation by h bits 

where 

h = offset + hits(m) + z-e (7) 

where offset is the bias offset (that is, the difference between biased exponent e and unbiased 
exponent e' in the general case where e^O), hits(m) is the number of bits in mantissa m, 
and z is the predetermined number of zeros added to binary number ms (step 506). In the 
implementation of Table 1, offset = 127, hits(m) = 23 , and z = 7 , yielding /? = 157 - e . The 

output of signed-shift-right shifter 1 10 is [xj in two's complement format. A signed-shift- 

right operation shifts the contents of the shifter by h bits while filling the new MSBs created 
by the shift with the MSB of the contents before the shift, as is well-known in the relevant 
art. In the implementation of FIG. 1, the MSB before shifting is the sign bit s. In one 
implementation, the signed-shift-right is performed by a signed-shift-right shifter 110. 

FIG. 6 is a logical block diagram of a format converter 600 according to an 
implementation of the present invention. FIG. 7 is a flowchart depicting a process 700 
executed by format converter 600 to convert floor ms fi-om conditionally-decremented 
floating-point format to floating-point format according to one implementation. 

Process 700 increments binary number by one if ^ = 1 (indicating that the floating- 
point number is negative), and does nothing otherwise, thereby producing an incremented 
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value (step 702). In the implementation of FIG. 6, the decrement operation is performed by a 
incrementer 602. In another implementation, as shown in FIG. 6, process 700 instead 
increments binary number ms by one if a conditional sign bit 55 = 1, and does nothing 
otherwise, thereby producing the incremented value (step 702). 

Process 700 then replaces the most-significant bit (MSB) of the incremented value 
with the exponent bits e and the sign bit 5, such that the sign bit s is the MSB, thereby 
producing the floor of the binary floating-point number in floating-point format. Li another 
implementation, as shown in FIG. 6, process 700 replaces the MSB of the incremented value 
with conditional exponent bits 63 and the sign bit 5, such that the sign bit s is the MSB, 
thereby producing the floor of the binary floating-point number in floating-point format. 

The present invention can be implemented in software. One implementation, in C 
code, is listed below. The code retums the floating-point floor in two's complement format 
as m6, retums the floating-point floor in floating-point format as B, and retums the fractional 
remainder as f9. 



// floor. cpp 
// 

#include "stdafx.h" 
#include <stdlib.h> 
#include <stdio.h> 
#include <math.h> // floor 

#define LIM(n) ( (n) >31 ? 31 : (n) <0 ? 0 : (n) ) 

ftdefine BlT(x,n) ( ( (x) >> (n) ) & 1) // x[n] 

#define BITS(x,m,l) ( ( (x) » (1)) & ones » (31- (m) + (1) ) ) // x[m:l] 
#define CC(n,x,Y) ( (x) « (n) | BITS (y, (n) -1 , 0) ) // ConCatinate 

{x,y} 

int main(int argc, char* argv[]) 

{ 

unsigned int fm, i, s, S3,x,m0,ml,m2,m3,m4,m7,m8,m9,ones=- 
l ; n , f m , s , s 3 , x , mO , ml , m2 , m3 , m4 , m7 , m8 , m9 , ones = - 1 ; 
int e^e3,m5,m6; 
float f ,f6,f7,f8,f9,fl0; 
int 

test [10] ={0x4f 000000, 0xb03 00000, 0x3f 000000, Oxbf 000000, 0, 0x800 00000, 0x80000 
001,1, 

Oxc0100000,0xc0000000} ; 
// for (i=0;i<10;i++) { 
// x = test [i] ; 

for (n = 0;n==0;n+=0x80000000) 

for (x=n;x<n+0x4f 000000, •X++) { f , f 3 , f 6 , f 7 , f 8 , f 9 , f 10 ; 
for (n = 0;n==0;n+=0x80000000) // test positive and neg sign 

for (x=n;x<n+0x4f 000000 ;x++) { // 
all values except NaN's 
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// FLOOR and REMAINDER begin HERE 



// Fig 1 Remainder 



f = * (float *)&x; 

mO = *(int *)&f; 

s - BIT(x,31) ; 

e = BITS(x,30,23) ; 

ml = CC(23, (e>0) ,BITS{x,22,0) ) ; 

S3 = s Sc (e>126) ; 

m2 = (ml - S3) ScOxf fffff ; 

fm = Oxffffff >> LIM(e-126); 

e3 = e<127 ? 127*s : e; 

m3 = m2 & -fm | (-s3) & fm; 

m7 = CC(31,s,CC(23,e,BITS(m2,22,0) ) ) ; 

m8 = CC(31,s,CC(23,e3,BITS(m3,22,0) )) ; 

f7 = * (float *)&m7; 

f8 = * (float *)&m8; 
* (float *)5cm9 = f7 - fS; 

f9 = * (float *)&m9; 



// FIG 4, INTEGER FLOOR 

m4 = m3 ^ (s * Oxffffff) ; 

m5 = CC (31, s,m4<<7) ; 

m6 = m5 >> LIM(157-e) ; 

// Fig 6, FLOATING POINT FLOOR 

*(int *)&f3 = CC(31,s,CC(23,e3,BITS(m3+s3,22,0) ) ) ; 

// TEST RESULTS AGAINST FLOOR IN C LIBRARY 
f6 = floor (f) ; 
flO = f-f6; 

if ((f9l=fl0) I (m6!-(int)f6)) (e < 158) & (m6 1 = (int) f 6) | (f3l=f6) ) 
printf ("\nx-%08x, ",x) ; // floor{-0.0) returns -1 

instead of 0 . 

if ( (x&0x07fffff)==0) 

printf (" . ") ; 

}} 

return 0; 

} 



The invention can be implemented in digital electronic circuitry, or in computer 
hardware, firmware, software, or in combinations of them. Apparatus of the invention can be 
implemented in a computer program product tangibly embodied in a machine-readable 
storage device for execution by a programmable processor; and method steps of the invention 
can be performed by a programmable processor executing a program of instructions to 
perform fimctions of the invention by operating on input data and generating output. The 
invention can be implemented advantageously in one or more computer programs that are 
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executable on a programmable system including at least one programmable processor 
coupled to receive data and instructions from, and to transmit data and instructions to, a data 
storage system, at least one input device, and at least one output device. Each computer 
program can be implemented in a high-level procedural or object-oriented programming 
5 language, or in assembly or machine language if desired; and in any case, the language can 
be a compiled or interpreted language. Suitable processors include, by way of example, both 
general and special purpose microprocessors. Generally, a processor will receive instructions 
and data from a read-only memory and/or a random access memory. Generally, a computer 
will include one or more mass storage devices for storing data files; such devices include 
10 magnetic disks, such as intemal hard disks and removable disks; magneto-optical disks; and 
optical disks. Storage devices suitable for tangibly embodying computer program 
instructions and data include all forms of non- volatile memory, including by way of example 
semiconductor memory devices, such as EPROM, EEPROM, and flash memory devices; 
ril magnetic disks such as intemal hard disks and removable disks; magneto-optical disks; and 

^',"15 CD-ROM disks. Any of the foregoing can be supplemented by, or incorporated in, ASICs 
ly (application-specific integrated circuits). 

^ ^ A number of embodiments of the invention have been described. Nevertheless, it will 

y be understood that various modifications may be made without departing from the spirit and 

C3 scope of the invention. For example, apparatus 100 and process 200 can easily be modified to 

find only one of the floating-point floor and fractional remainder, rather than both. Further; 

apparatus 100 and process 200 can easily be modified to handle numbers having a different 

number of bits than those described above. Accordingly, other embodiments are within the 

scope of the following claims. 
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